BINOMIAL FORMULA FOR MACDONALD POLYNOMIALS 

AND ITS APPLICATIONS 
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Andrei Okounkov 

Abstract. We generalize the binomial formula for Jack polynomials proved in 
[002] and consider some applications. 



§1 Binomial formula 

Binomial type theorems (that is Taylor and Newton interpolation expansions 
about various points) are powerful tools for handling special functions. The high- 
school binomial formula is the Taylor expansion of the function fix) = x l about 
the point x = 1. Its (/-deformation is the Newton interpolation with knots 



x = l,q,q , ... 



which reads 
(1.1) 

where 



x 



/ 

m 



(x-l)...(x-q m - 1 ), 



I 

m 



[q l -1) ...{q l -q™- 1 ) 
-,m—l s 



J(/ (<r - 1) . . . (<r - ^ . 

is the (/-binomial coefficient. Denote the Newton interpolation polynomials in the 
RHS of (1.1) by 



Pt{x-q) 



fx -11 



[x — q 



fc = 0,1,2, 



The formula inverse to (1.1) is the Taylor expansion of Pj*(x;q) about the point 
x = 



;i.2) 



(x-l) 



[x — q 



Ml-l)/2 



B-i: 



l—m 



I 

m 



x 



m(m-l)/2 



which is essentially the fundamental (/-binomial theorem [An,GR]. 



The paper was written during the author's stay at the Institute for Advanced Study in Prince- 
ton. I am grateful to IAS for hospitality and to NSF for financial support (grant DMS 9304580). 
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The formulas (1.1) and (1.2) are limit cases of the formula 



(1.3) 



Pj*(ax; q) 



a 



I 

m 



P*(x;l/g) 

l/q P ^Q) 



as a — > oo and a — > respectively. Since 



/ 

m 



= p* (<?-'; i/g) 

i/, P*(<r m ;l/«) 



the formula (1.3) is symmetric in g z and x. The formula (1.3) is easy to remember 
in the form 



(1.4) 



l« = j:(...)^( ! -';l/r B (x;l/ g ) ) 



where the factor 



( ■ ' ■ ) p* 



P£(g-™;l/g)P£(a;5) 

can be easily reconstructed using that 

(1) this factor does not depend neither on x nor on /, and 

(2) the highest degree term of Pj*(x; q) is x . 

In this paper we consider a multivariate generalization of the formula (1.4) for 
symmetric functions. Let q and t be two parameters. Let fx be a partition. Recall 
the following notations of Macdonald [M] . Set 



n(/x) = ^(z-l)^ = ^^(^.-l)/2. 



For each square s = £ /U the numbers 



a(s) = m-j, a'(s)=j-l, 
l(s) = (jfj — i, l'(s) =i — l, 

are called the arm-length, arm-colength, leg-length, and leg-colength respectively. 
We shall need the following normalization constant 

(1.5) H(n; q, t) = r 2n ^q n ^ J] (y^H Y (s) - l) . 

sen 

Denote by V the set of all partitions. This set is partially ordered by inclusion 
of partitions fx C A. Consider the following multivariate Newton interpolation 
polynomials 
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Definition ([Kn,S2,Ok3]). Suppose fx has < n parts. The interpolation Macdonald 
polynomial 

P^ (x\, . . . , x n ; q, t) 

is the unique polynomial in x\, . . . , x n such that 

(1) P*(xi, . . . , x n ; q, t) is symmetric in variables x$~ l , 

(2) degP; = \fil 

(3) P;(q»;q,t) = H(fi;q,t), and 

(4) P*(q x ; q, t) = for all partitions A such that fx <f_ A. 

Here P*(q x ; q, t) = P*(q Xl , . . . , q Xn ;q, t). Another name for the same polynomials 
is the shifted Macdonald polynomials. 

The connection with ordinary Macdonald polynomials will be explained below. 
It is easy to see that the polynomial P*(x; q, t) is indeed unique provided it exists 
(which is not evident). It follows from the definition that P*(x;q,t) is stable, that 
is 

P*{x x , . . . , x n , 1; q,t) = P*(xi, . . . , x n ; q, t) 

provided fx has at most n parts. Using this stability one can define polynomials 
P*(x;q,t) in infinitely many variables. 

Remark. The condition (4) in the definition of P*(x;q,t) can be replaced by a 
weaker condition 

(4') P*(q x ; q, t) = for all partitions A ^ fx such that |A| < \fx\ . 

The first proof of the fact that (4') implies the property (4), which is in that case 
referred to as the extra vanishing, was given in [Kn]. 

The interpolation polynomials in the Schur functions case were studied in [00] 
(and, from a different point of view, in earlier papers cited in [00]). They are called 
shifted Schur functions, have many remarkable properties, and can be interpreted 
as distinguished basis elements of the center of the universal enveloping algebra 
W(gl(n)) (see also [Okl,2] where explicit formulas for the corresponding central 
elements were given). We have also conjectured (unpublished) the general combi- 
natorial formula for the polynomials P*(x;q,t) (the formula (1.6) below). Then 
F. Knop and S. Sahi proved a series of results for general (q, t) in [KS,Kn,S2]. In 
particular, it was proved that the polynomial P^(x; q, t) exists and has the funda- 
mental property (1.8) below. They used difference equations for the polynomials 
P*(x;q,t). Later the combinatorial formula (1.6) and also a (/-integral representa- 
tion for P^(x; q, t) were proved in [Ok3]. 

We call a tableau T on fx a, reverse tableau if its entries strictly decrease down the 
columns and weakly decrease in the rows. By T(s) denote the entry in the square 
s E fx. We have [Ok3] 

(1.6) P;{x;q,t) = J>tM) Ht lT ^ (x t(s) - /<->*-''<•>) , 

where the sum is over all reverse tableau on \i with entries in N and ifiriq, t) is the 
same (q, t)-weight of a tableau which enters the combinatorial formula for ordinary 
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Macdonald polynomials (see [M],§VI.7) 

(1.7) P^(x;q,t) = ^2^ T (q,t)Y[x T{s) . 

The coefficients ^(g, t) are rational functions of q and t. Explicit formulas for 
them are given in [M], (VL7.11/). It is clear from (1.6) and (1.7) that 

(1.8) P^{xi, ■ ■ ■ , x n ; q, t) = Pn(xi, X2t _1 , . . . , x n t 1 ~ n ) + lower degree terms . 

In other words, ordinary Macdonald polynomials are the top homogeneous layer of 
the interpolation Macdonald polynomials. 

It follows from the combinatorial formula (1.6) that 

(1.9) P;(o, • • • , a; q, t) = P M (1, t~\ . . . , t 1 ^- q, t)]j(a- qr»'«f-''W) , 

sen 

where a is repeated n times. Recall that (see [M], (VI.6.11')) 

i „a'(s)j.n— l'(s) 
, 1 i_ n s _ .nM + lftld-n) TT 1 ~ g ^ 

Now we can state our main formula which is the following generalization of (1.4) 
Binomial Theorem. 

P*(a,...,a;q,t) ' ^ P;(g-M; l/t) p*( a , . . . , a; g, *) 

Remarks. 

(1) By definition of P*(x;q,t) only summands with \x C A actually enter the 
sum in the binomial theorem. 

(2) By (1.5) and (1.9) all denominators in the binomial theorem have explicit 
multiplicative expression. 

(3) One can easily remember the binomial formula in the same way as the 
formula (1.4). 

(4) The binomial formula calls for notation 



qjt P*(q»;q,t) 



One can rewrite the binomial formula as an identity for the above (q, t)- 
binomial coefficients; see the formula (2.5) in the next section. Note that 
unlike (/-binomial coefficients the above (q, t)-binomial coefficients are not 
polynomials 1 . 



1 M. Lassalle recently pointed out that in particular case when X/^i is a single square these 
binomial coefficients appeared in [Ka2]. 
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(5) By (1.9) we have 

P*{a, . . . , a; q, t) = (-a) W^V^-^-^P^l/a, . . . , 1/a; 1/q, 1/t) 



(1.10) 



^ V ^,, . . . , ^j, „j y ^ ^ ^ ^ , 

Therefore we can rewrite the binomial formula as follows 
P x *(ax 1 ,...,ax n ;q,t) _ 



P*(a,...,a;q,t) 

y ( , M I t^) P;(Q- A ;l/g,l/t) P;(x n ,...,x i; l/g,l/t) 
^ ' q»W P;(g-M; 1/q, 1/t) P;(l/a, . . . , 1/a; 1/q, 1/t) ' 

Using (1.8) and the well known equality 

P x (x;q,t) = P x (x;l/q,l/t) 

we obtain as limit cases as a — > oo or a — > binomial formulas involving ordinary 
Macdonald polynomials. We have 

( , Pxjxi,..., x n t x ~ n - q, t) _ P;{q X ; q, t) P* (x u . . . , x n ; q, t) 



P A (1, . . . , t 1 — ; q, t) ^ P;(q»- q, t) P M (1, . . . , t 1 — ; Q, t) 

A 1 



and 



P*(q;i, ■ ■ ■ , x n ; q, t) = P;(<T A ; 1/g, 1/t) P^ . . . , g, *) 

l ' } P A *(0,...,0;Q,t) ^P;(g-M;l/g,l/t) P;(0,...,0;Q,t) ' 

as a — > oo and a — ► respectively. 

If we set t = q d and let q — > 1 then (1.11) turns into the binomial formula for 
Jack polynomials proved in [002]. Thus, the above binomial theorem is a (far 
reaching) generalization of that formula. Various results on binomial coefficients 
for Jack polynomials were obtained in [Bi,FK,Ka,Lasc,La,00]. 

In the next section we list some applications of the binomial formula. Those 
applications are quite different from the applications of the Jack degeneration of 
the binomial formula considered in [002,004,Ok5,Ok6] (see also [KOO]). Then in 
section 3 we recall some results of F. Knop and S. Sahi. The proof of the theorem 
is given in section 4. Finally, in Appendix, we discuss an elementary algorithm for 
multivariate symmetric Newton interpolation (to avoid possible confusion: this is 
not the interpolation considered by Lascoux and Schiitzenberger.). 

A binomial formula for the Koornwinder polynomials (which generalize Macdon- 
ald polynomials for the classical root systems) was proved in [Ok4]. 

I would like to thank G. Olshanski and S. Sahi for helpful discussions. 

§2 Applications 

Consider the following matrix whose rows and columns are indexed by partitions 
with at most n parts 



S(a,n;q,t) = 



P^aq-^,...,aq-^;q,t) 
P*(a,...,a;q,t) 



From the binomial theorem we have 
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Symmetry. The matrix S(o, n; g, t) is symmetric for any a. 

This symmetry is a simple yet central result. Since it interchanges the label and 
the argument of the polynomial P* all properties of P* come in pairs, where one is 
obtained from the other using the above symmetry. For example, one checks that 
we have two following pairs 



combinatorial 
formula 



g-integral 
representation 



difference 
equations 



Pieri-type 
rules 



(2.1) 



Similarly, it is evident from (1.11) that the matrix 

P x (q-\...,q^t 1 - n ;q,t) 



S(oo,n;l/g,l/£) 



P x (l,...,t^;q,t) 



is symmetric, which was conjectured by I. Macdonald, proved by T. Koornwinder 
(unpublished, see [M], section VI. 6), and, by different methods, in [EK] and [Ch]. 
The binomial formula is just the Gauss decomposition of the symmetric matrix 

(2.2) S(a, n; q, t) = B'(l/g, 1/t) A(a, n; 1/q, 1/t) B(l/g, 1/t) 

where prime stands for transposition, the matrix 



(2.3) 



B(q,t) = 



Q,ty 



is upper triangular (with respect to the partial ordering of partitions by inclusion) 
and unipotent, and the matrix 



A(a,n;q,t) = ( 6 X »(-1) 



q' 



Km') 



t n M p*(l/ a ,...,l/a;q,t) 



At, A 



is diagonal. 

The symmetry of the matrix S has the following combinatorial interpretation. 
Suppose that for some b E N 

A, v C (3 

where (3 is a box f3 = (b n ). Denote by 

p\X=(b-X n ,...,b-X 1 ) 
the partition complementary to A in (3. Setting a = q b we obtain 

PZ(qK";q,t) _ P* v {q p \\ q,t) 



(2.4) 



P* x (qP;q,t) Pi(q^,q,t) 
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The skew diagram ((3\u)/X is depicted in the following Fig. 1 (where the diagram 
v is rotated by 180°) 



A 



v 



P 



Fig. 1 

Introduce the following trinomial coefficients 



" p ' 




'P 






A, v 


q,t 


V 


q,t 


A 



q,t 



Using them one rewrites (2.4) as follows 
Symmetry II. 



In particular, 



" p ' 




' P " 


A, v 


q,t 


v, A 



q,t 



~p~ 




P 




V 


q,t 


_p\v_ 


q,t 



In the shifted Schur function case the symmetry is obvious (see Fig. 1) from the 
following formula (see [00], §8, and also [002], §5 for Jack generalization) 



dim(/?V)/A, 



where dim stands for the number of standard tableaux on a (skew) diagram. 
Observe that 



" p ' 




f \P\ N 


\ dim A dimz/ 


A, v 


1,1 




' dim/3 



A{q\n-MqMt)= U v (-1)^ ! 



~p 


1 ) 




l/q,l/tj 



Therefore, the binomial formula can be rewritten as follows 
Binomial Theorem II. 



(2.5) 



' p ' 




'P 




'P 


A, v 


q,t 


A 


q,t 


V 



q,t 



Et- 1 ) 1 



7 n(/x') 



A 




V 




P 


— 1 








l/q,l/t 




l/q,l/t 



Note that only summands with fx C A fl v actually enter this sum. 
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We denote by A*(n) the algebra of polynomials in x\, . . . , x n which are symmetric 
in variables 

We call such polynomials shifted symmetric. The polynomials P*(x;q,t), where fx 
ranges over partitions with at most n parts, form a linear basis in this algebra. 

The binomial formula gives transition coefficients between various mutually tri- 
angular linear bases of A*(n). Combining, for example, (1.11) and (1.12) and taking 
into account that by (1.9) 

P*(0,...,0;q,t) = ( )M n ( M ')+-n(M) 
P M (M-i,...,ti-»; ? ,f) 1 ] q 



we obtain 



+n(v) Q n(p') 

^ = (-i) M ^£(-i) l "'W 



A" 




n 






q,t 


V 


1/9,1/* 



This is equivalent to the following proposition 

Inversion. The inverse matrix of the matrix (2.3) is given by 

(2.6) B- 1 ^,*) = C- 1 ^,*) B(l/q, l/t) C(q,t) , 

where C(q,t) is the diagonal matrix 

n(>') 



C(q,t)= 18^(-1)^ (J 



/i, A 



Note that the inverse of the (truncated) matrix (2.1) in the case when q is a root 
of unity was calculated by A. A. Kirillov, Jr. [Ki] and also by I. Cherednik [Ch]. 

The algebra A*(n) is filtered by degree of polynomials. We denote by A* the 
projective limit of algebras A*(n), n G N with respect to homomorphisms 

f(x u . . . , x n+l ) i-> f(x 1 , ...,x n ,l). 

Note that f(q x ) is well defined for all / G A* and A G V. Consider the Newton 
interpolation of a polynomial / G A* with knots 

x = q X , XeV, 

or, in other words, the expansion of / in the linear basis {P*{x; q, t)}^ e <p of A*. 
The coefficients /(//; q) of this expansion 

f(x) = J2tt^q) p ^^t) 

can be found from the following non-degenerate triangular (with respect to ordering 
of V by inclusion) system of linear equations 

(2-7) f(q X ) = J2f^ti P ^^ t ^ 

A*CA 

Using (2.6) we can solve the equations (2.7) and obtain 
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Explicit Newton Interpolation. The coefficients f(n;q) are given by 



^ 1 / 1 n (/ i l v ) 

< 2 ' 8 > /(«9) = 7W^T7TE(- 1 )"' /H l 



P*(q»:q,t) ^ ' t n ^/ v ) 



v 



l/q,l/t 



where n(fx/u) = n(fx) — n(u). 



In particular, this formula can be used for expansion of any homogeneous sym- 
metric polynomial in ordinary Macdonald polynomials. 
One can think of the function 

f:fi^f(fi,q) ^eV 

is an of analog of the Fourier transform of the function 

/ : /* ~ /(«") ^ev. 

The inverse transform is obviously given by (2.7). 

Combining the formula (2.8) with the Knop-Sahi difference equation for the 
polynomials P* one can obtain a nice multivariate generalization of the classical 
algorithm for Newton interpolation with knots 0, 1, 2, . . . (see Appendix). 

The binomial theorem provides a natural proof of the (/-integral representa- 
tion for interpolation Macdonald polynomials found in [Ok3] . In [Ok3] the integral 
representation was guessed and checked. One can deduce the integral representa- 
tion from the formulas (1.7) and (1.11) in exactly the same way as it was done in 
[002] for Jack polynomials. 

Notice that the the binomial formula makes certain sense for A which is not a 
partition. An example of such analytic continuation is the formula (2.10) below. 
Notice, however, that the RHS may or may not converge to the LHS in the usual 
analytical sense. 

As two particular cases of the binomial formula, one obtains two following gener- 
ation functions for the 1-column and 1-row interpolation Macdonald polynomials. 
By definition, set: 

e k( x i t) = P(ikj( x 'i q-, t) , 

h* k (x;q,t) = P* k) {x;q,t) , fc = l,2, 

The polynomials e k (x;t) do not depend on q for it follows from the combinatorial 
formula (1.6) that 



e* k (x;t)= t k -^^Y[( Xl -t s - k ) 



ii<---<ik 



Let u be a parameter. We have the following 
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Generating functions. 

(2-9) 

1 + Xit x ~ % ju _ ^ e l( x '^) 



(2.10) 



{ . + *!-*/« ^ (U+1)(U+ 1 /*)... 



(xi*-* /w;g)oo (* 1_i / u ;9)oo ^ (g;g)fc (w - l)(w - g) . . . (w - g^" 1 ) ' 

Here (z; q) r = (1 — z)(l — qz) ... (1 — q r ~ 1 z). Note that above expansions make 
sense as formal power series in 1/u with coefficients in A*. 

This generating function (2.9) is a minor modification of the generating function 
in [00], theorem 12.1. The analog of the generating function (2.10) for shifted Jack 
polynomials was found by G. Olshanski (unpublished) and it was used by him to 
prove the Jack degeneration of the 1-row case of the combinatorial formula (1.6). 

It is clear that if 

X n+1 = X n+2 = • • • = 1 

then (2.9) is just the binomial formula for A = (l n ) and a = —t 1 ~ n /u, whereas 
(2.10) is the binomial expansion of 

PZ(x n /u, . . . , xi/u; 1/q, l/i) 
P*(l/u,...,l/u;l/q,l/t) 

analytically continued to the point A satisfying 

q X = (l/t,...,l/t). 

It is easy to check both expansions directly. 2 Since it will be convenient for us to 
use (2.9) below in section 3, we shall give such a direct proof, which is essentially 
the argument from [OO], section §12. Same argument applies to (2.10) 

Proof of (2.9). Let e° k {x;t) be the coefficients in the expansion 

(2U] TT l + Xjt^/u _ y> e° k (x;t) 

K ' ' l + t^/u f^ Q (u + l)(u + l/t)...(u + l/t k - 1 )' 

We want to prove that e k (x;t) = e k (x;t). It is clear that e° G A*. Set Xi = ^uz, 
where and z are new variables and let u — > oo. One obtains 

e1(x; t) = ek(xi, xit~ X . . . ) + lower degree terms . 



2 In fact, one can generalize the expansion (2.9) by replacing the sequence 1, l/i, l/i 2 , ... by an 
arbitrary sequence. That type of expansions was considered by A. Abderrezzak in [Ab]. However, 
only very few facts about interpolation Macdonald polynomials admit such a vast generalization, 
see [Ok7]. 
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Therefore, it suffices to prove that 



e° n (x u . 

Substitute x n = x n +i = 



,x n _i,l,...,l;t) = 0, re=l,2,.... 
• = 1 in (2.11). Then the RHS becomes 



n-l 

n 

i=i 



U + xd 



1-i 



U + t 



1-i 













q,t 


y. 


l/t,l/q 



In its expansion of the form (2.11) all summands with k > n vanish. □ 

The expansions (2.9) and (2.10) are related by the duality for the interpolation 
Macdonald polynomials, see [Ok3], theorem IV. That duality (in particular, the 
proposition 6.1 from [Ok3]) implies the following duality for (q, £)-binomial coeffi- 
cients 



(2.12) 



§3 Difference equations for P*(x;q,t). 

In this section we recall some results due to F. Knop and S. Sahi [Kn,S2]. Only 
the explicit formulas for higher order difference operators are new. 
Denote by T XjQ the shift operator 

[W] (*) = /(«*). 

Observe that for any polynomials pi ,p% in x such that 

Pl(0)=P2(0) 

the operator 



Pi(x)f(x) -p 2 (x)f(qx) 

x 

maps polynomials to polynomials. Set 

Ti = T Xu i/ q . 

Consider the following n-dimensional difference operator 



D{u) = V(x) _1 det 



{1 - UXl \l-t Xl )^- {1 - Xl)3 T, 



Xi 



Here u is a parameter and V(x) is the Vandermonde determinant. The determinant 
is well defined since multiplication by Xi and the shift Tj commute provided i ^ j. 
It is clear that D(u) maps A(n) to A(n) and does not increase degree. 
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For any subset Id {1, . . . , n} set 

iei 

Using the Vandermonde determinant formula one expands D(u) as follows 

D(u)=i[x- 1 £ (-lyk-iki-xjiki-vxi) n x -0^ t ^ 

/C{l,...,n} iS/ ig/ i€l,j<£l % 3 

where the omitted exponent of t equals (n — \I\)(n— \I\ — l)/2. 
Now shift the variables 

and consider the operator 

D*(u) = Hx- 1 J2 t^-W^-W- 1 )' 2 

JC{1,...,»} 

nd -x.r-i n<i ii T^^ ' 

ie/ igj iei,j<£i 1 3 

which maps A*(n) to A*(n). 

Below we shall have to evaluate 

[D*(u)f](x), feA\n) 

at the points of the form 

(3.1) x={aq^,...,aq^), &eZ, 

where £i > • • • > £ n an d a is a parameter (which will be the same parameter a as 
in the binomial theorem). Fix a subset I and set 



& - 1, i e * 
6, H 1 - 



We have the following obvious 
Lemma. Let x be as in (3.1). Then 



n j 

- L - L Xi — Xnt 3 

3 



i#ei>--->e; 

This lemma results in the following 
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Theorem. 

(3.2) D*(u)P;(x;q,t)= ^[(q-^f' 1 - ut^ 1 ) j P;(x;q,t). 

This theorem (without explicit formula for D*(u)) is due to F. Knop and S. Sahi 
[Kn,S2]. 

Idea of Proof. Using that 

JJ(1 - x % t n ~ % ) =0 if x n = 1 and % E I 

iei 

and the lemma one easily checks that 

[D*(u)P;](q x ) = 0, AeP,A^ 

and that 

[D*(u)P;] (<f) = ^[(q-^t* - 1 - ut"- 1 )^ P;{q») . □ 

In particular, the operators D*(u) for different values of u commute. Consider 
the algebra 3) generated by these commuting difference operators. 

This commutative algebra is isomorphic to the algebra A 1 / t (n) under the Harish- 
Chandra map which maps an difference operator D e 3) to the polynomial d G 
A 1 /*^) such that 

DP;(x;q,t) = d(q-»)P;(x;q,t). 
Let .Dfc be the coefficient 

D*(u) = ^(-l)"-*^- 1 )^-*)^ - - ...(u- l/t n - k - 1 ) D k 

k<n 

in the Newton interpolation of D*(u) with knots 

u = 1,1/t, 1/t 2 , ... . 

By (3.2) and (2.9) we have 

D k P; (x; q, t) = e* k (q^; 1/t) P; (x; q, t) . 
The operators D k are generators of the algebra D. 
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§4 Proof of the binomial theorem. 

We begin with general remark about the proof. As it is usual in the Macdonald 
polynomials theory, various pieces of the puzzle can be arranged in various combi- 
nation to give a proof. Below we shall use the approach based on the Knop-Sahi 
difference equations. However, it will be clear that the formula (4.4), which states 
that the LHS is a polynomial in q~ x with a certain type of symmetry, is the key 
point of the proof. Instead of using difference equations, one can deduce (4.4) 
and hence the binomial theorem from the g-integral representation of interpolation 
Macdonald polynomials [Ok3]. 

In this section by c 7 we always denote some coefficients. This notation is used 
to show which summands possibly enter the formula. 

It is clear from definitions that 

p*=E(-») m «-n (1 ~f ~' } ii x, ~-?#-? T ' + 

\I\=k i€l % i€l,J0 1 3 \I\<k 

Suppose 

x=(aq^,...,aq^), £eZ, 
where £i > • • • > £ n . Then by the lemma 

(4.1) [D k f\(x) = '£ t c ri f(aq r ",...,aq r ">) 

v 

with some coefficients c v , where 

Vi > • • • > Vn , 

(4.2) £_^ e {o,l}, i = l,..., n , 

Order the n-tuples r] lexicographically as follows 

7] > fj 

if rj n < fj n , or if r] n = fj n and i] n -i < rj n -i, and so on. Then the maximal 77 
satisfying (4.2) is 

Vmax (£l? • • • ? £,n—ki Cn—k+1 !>•••> £n 1) • 

It is clear that the coefficient 

in (4.1) is a non-zero rational function of a. 

Now let fx be a partition and let fi' be the conjugate partition. Iterating the 
above argument we see that 



(4.3) 



U D »f 



{a,...,a) = ^c v f{aq Un ,...,aq Vl ) , 
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where c v are some coefficients and v ranges over partitions such that \v\ < and 

V<fl' 

in the lexicographic order of partitions. Observe that the coefficient in (4.3) 
does not vanish. Therefore there exist an operator D(ji, a) e D 

D(n,a)= c »n D v'i 

\u\<\n\,v<tJ, i 

such that 

[£>(//, a)/] (a,...,a) = f(aq-^,...,aq-^). 
Denote by e/ M)Cl G A 1 / t (n) the Harish-Chandra image of D(n, a) 

D{n, a) P{{x- q, t) = d^ a {q~ x ) P* x {x- q, t) . 

By definition of D(fi, a) we have 

Since the degree of e£ is we have 

degrf M , a < \n\ . 
Now consider the Newton interpolation of the function 

= P^xi,...,ax n ;qt) A t (n) 

M J PJ(a,...,a;g,*) 1 J 

with knots 

This interpolation has the form 

(4-5) = E W • • • , X i; 1/t) . 

The coefficients fr^A, a are linear combinations of the values 

PJ(ag- / »,...,ag- ,/1 ;?,*) 



P£(a,...,a; 

By (4.4) we have 



& M ,A,a = ^,a(? A ) 
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for certain polynomials 

& M , a G A 1/4 (n) , deg& M>0 < |//|. 
Next observe that the highest degree term of the LHS of (4.5) equals 

IAI 



a 



P*(a,...,a;q,t) 



P\(x 1 ,x 2 t 1 ,...,x n t 1 n ;q,t) 



and 



P;(x n , . . . , xi; 1/q, 1/t) = t^-VMp^xu x 2 t~\ x n t x ~ n - q, t) + . . . , 
where dots stand for lower degree terms. Therefore 

, |A| < and fx ^ A , 



a lMlt(i-n)lMlp*( a; . . . , a; t)' 1 , (j, = A . 



Since degfr^a < the polynomial is uniquely determined by its values at the 
points x = q~ x , |A| < This implies 



aMt (i-n) M p*{x-l/q,\/t) 



P*(a,...,a;q,t)P£(q-»;l/q,l/t) 
and concludes the proof of the theorem. 

Appendix. Symmetric Newton's algorithm 

The formula (2.8) can be turned into a quite efficient algorithm for multivariate 
symmetric Newton interpolation. 

Reccurence relations for generalized binomial coefficients were first used in the 
Jack polynomial case by M. Lassalle, see [La] 3 . Those relations correspond to a 
Pieri type formula for interpolation Macdonald polynomials (see [002], §5). They 
are related to the Knop-Sahi type relations used below by the symmetry, see §2. 

Let A be a diagram and suppose that v is obtained from A by removing a corner. 
We shall denote this by writing 

v/\. 

Denote the corner being removed by see Fig. 2 



v 



r 



X 



Fig. 2. The diagram v is obtained from A by removing the corner 



3 See also M. Lassalle's recent preprint Coefficients binomiaux generalises et polynomes de 
Macdonald, S10. 
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Consider the outer corners of the diagram v (denoted by o in Fig. 2) and also the 
inner corner of v excluding the one where the square \jv is attached (those inner 
corners are denoted by • in Fig. 2). By definition, set 



inner corners 



Jl (l - g Ay(*,o) t Ax(o,*) j ' 



outer corners o 



where 

Ay (it, •) = y-coordinate ( * ) — y-coordinate ( • ) , 
Ax(», it) = x-coordinate ( • ) — x-coordinate (it) • 

Now let / G A* be the shifted symmetric polynomial to be interpolated. For any 
pair n C A set by definition 



F(/,,A) = (-i) 



n(A'/ M ') 



fn{X/n) 



l/q,l/t 



where n(\/(i) = n(X) —n(fi). Then, the interpolation formula (2.8) can be rewritten 
as follows 



(5.1) 



/(a) = £fgm) 



On easily computes that the Knop-Sahi difference equation 



implies that 



(5.2) ^ W (z//A; ? ,t) 



V 




A 








A 




l/q,l/t V 


1 


q,t 


1 







l/q,l/t 



Therefore, the function F(ji, A) satisfies the following reccurence relation 



F(n,\) = - ^ i— ^W)> 



A" 




[I 


1 


q,t 


1 







7l'(A/i/) 

which, together with the obvious initial condition 

*W) = /(g"), 
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more simple than the coefficients 



give a simple algorithm for calculation of this function. The summation (5.1) can 
be performed simultaneously with the computation of F(fi, A). 

Both in formula (2.8) and in the above algorithm one has to run a cycle over 
all pairs of partitions n C A. This means of order N 2 iterations to interpolate in 
N points. (The same number of iterations is sufficient for direct solution of the 
equations (2.7).) However, the above algorithm has one clear advantage: 

The function u(v/*\;q,t), which we have to compute on each step, is much 

.In other words, one can expand a 

polynomial / in interpolation Macdonald polynomials without actually knowing 
those quite complicated polynomials explicitly. 

In the Jack limit this algorithm can be made exponentially faster. Using the fact 
that the LHS of (5.2) depends in the Jack limit only on \/i\, not on particular form 
of /i, one is able to use on each step some previously computed sums as follows. 

We shall use the parameter 9 as in [002]. Let 

f(x 1 ,x 2 ,...) 

be a polynomial, symmetric in variables Xi — 6i. By definition, set 



The polynomials 



"A" 




lim 


"A" 






9 




.P. 






X 








A 


9 



form a linear basis in the algebra of polynomials symmetric in variables Xi — 6i. 
The coefficients /(//) in the expansion 



/ = £/(/*) 



can be found from the following specialization of the formula (2.8) 

A 



/(A) = £(-l) |A/ "' 

A*CA 

Now, instead of computing the numbers 



/(/*)• 



e 



F(^,A) = (-1)I a /mI 
we shall deal with the sums 
FW(A)= £ 



MCA 

\\/»\<k 



k 

|A//xj 



MCA 
|A//i|<fc 



k 

|A/^| 



/(/*), 
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where 

fe = 0,...,|A|. 

Then we have 

F<°>(A) = /(A), 
^ (|A|) (A) = /(A), 

and the following reccurence relation 

ir(*+i)(A) = F«(A) - t^t V u{v/\\ 6) F^ k \u) , 

A 

v/X 

where (see Fig. 2) 

oj{v/\; 9) = lim u{v/\; q, q e ) 

J] (Ay(*,.) + 0Ax(.,*)) 

inner corners • 

~ II (Aj/(*,o) + 6IAx(o,*)) ' 

outer corners o 

The advantage of this algorithm is that instead of a cycle over all pairs of partitions 
C A one has a cycle over all partition A and all numbers k, such that < |A|, 

which is much faster. 

In the case of one variable this algorithm specializes to the classical Newton's 

algorithm for interpolation with knots 0, 1, 2, . . . because 

u{{l-l)/{l);9) = l. 
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